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ABSTRACT 

A flight loading profile based solely on accelerometer 
data Is developed which accounts for all fatigue damaging 
loads experienced by high performance aircraft. A random 
sequence of the flight loads Is generated In a computer 
program that applies the ESDU Cumulative Damage Hypothesis 
for damage prediction of fleet aircraft. Data from actual 
accelerometer readings are used In evaluating vjhether 
reasonable results can be obtained from the program. The 
feasibility of Implementing the program is discussed along 
with recommendations improving the suitability of the 
program for m.ore accurate predictions. 
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I. INTRODUCTION 



The design phase of new aircraft involves the use of 
an estimated flight loading profile in order to make a ser- 
vice fatigue life estimate of the aircraft. A full-scale 
fatigue test is normally performed on the aircraft in the 
early stages of production to verify the usable life predic- 
tion; however it is inevitable that the actual loads the 
aircraft experiences throughout its service life will be 
different than the design spectrum. The fatigue monitoring 
program presently used for naval aircraft attempts to record, 
without regard to sequence, the actual loads an aircraft 
experiences by means of a four level exceedance accelerometer. 
Analysis of flight data from oscillograph recordings indicates 
that each type of aircraft has a distinguishable flight load- 
ing profile and it has been determined that the order in 
which these loads occur on a particular type of aircraft 
has a significant effect on the damage predlciton. This 
thesis proposes a procedure vjhich estimates the actual flight 
profile from accelerometer data and statistically formulates 
a probable order in which these loads may have occurred. A 
fatigue cumulative damage theory is then incorporated v.l-iich 
accounts for the sequential effects of loads in fatj.gue damage. 
The method is adapted to the computer v/hich is necessai:'y for 
monitoring the fatigue damage of a .large fleet of aircraft 
and provides a damage calculation for each aircraft in 
approximately 2.5 seconds. 



B 



II. DAMAGE PREDICTION 



Structural components which experience repeated loadings 
are susceptible to fatigue damage for loading amplitudes 
v;hich occur well below the fracture strength of the material. 
These loading conditions occur frequently on aircraft, since 
the requirement of maintaining a particular flight pattern 
requires many loading changes. As a result, the design of 
aircraft structural components must account for the damaging 
effects of these repeated loadings, which necessitates an 
accurate means of predicting the fatigue damage accumulated 
due to these loading conditions . 

Fatigue life will be understood to be the life until a 
visible crack is present, or the life until complete failure 
of a small component where the crack propagation phase 
v/ould be short in comparison to the crack initiation period. 
Damage v/ill refer to the percent of fatigue life expended 
due to a prescribed loading condition. Methods for predic- 
ting damage invariably rely upon some application of the 
time-tested Palmgrem-Miner Rule, better known simply as 
Miner's Rule. 

A. MINER’S RULE 

Both Palmgren and Miner are credited vjith introducing 
what has become the most well known Incremental damage 
theory [Ref. 1]. If a component has experienced n^ load 
cycles v;ith the .same mean load and load ainjjliLude, tl:o 
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portion of the fatigue life expended is n^/N where N is the 
life to failure in a constant-amplitude test with the same 
mean and amplitude. Failure is assumed to occur if the sum 
of the consumed life portions equal 100 percent. This 
implies that the condition for failure is: 

En./N = 1 (1) 

1 

This failure criterion was extended to a component 
experiencing load cycles of various amplitudes and means 
such that the condition for failure is : 

E(n^/N^) = 1 (2) 

Specimen testing of different materials under varying 
mean and amplitude loads provide the to use in Miner's 
Rule. The family of curves developed from specimen testing 
for a particular stress concentration factor are referred 
to as S-N curves. 

1 . Advantages 

Due to the simplicity of application. Miner's rule 
has enjoyed v/ldespread popularity as a means of predicting 
damiage in aircraft. It is particularly suited to a large 
fleet of aircraft where loading data obtained from acceler- 
ometer readings can be readily stored in computers, which 
allov/s for fully automatic processing of data o)id evaluation 
oi’ damage prediction for individual ali-craft. Minei-'s rule 



10 



provides a quick "rule-of-thumb" indication of fatigue life 
expended since as the number of cycles of loading Increases 
a greater portion of fatigue life is expended. 

2 . Limitations 

Due to the complexity of aircraft loading cycles, 
Miner’s rule has proven to be a rather conservative estimate 
of fatigue damage. For Instance Lambert [Ref. 2] applied 
a typical aircraft block program to a rlvetted lapjoint of 
2024-T3 clad sheet aluminum. The ratio of actual life to 
Miner's estimated life obtained was 2.64. Through an exten- 
sive survey of the area of fatigue damage Schijve [Ref. 1] 
concludes that a Miner's Rule approximation provides a con- 
siderably inaccurate estimation of the fatigue life expended 
by aircraft. Economically it is not sound to retire an 
aircraft from service when actually only 37 percent of its 
fatigue life has been expended. 

B. SEQUENCE EFFECTS 

Since loading amplitudes vary randomly on an aircraft 
throughout its life, it has been determined that the order 
in which these loads are applied must be significant in 
the fatigue prediction. S-N data applied to Miner's Rule 
are based on data points obtained by specimen failure at a 
constant amplitude for a specific mean stress; therefore, 
it is highly unlikely that damage accumulation in aircraft 
components, vfhlch experience a wide range of amplitudes from 
cycle to cycle, v;ould accurately adhere to a Miner's Rule 
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approximation. The order of loads is not taken into account 
by Miner’s Rule. Schljve [Ref. 1] attributes the difference 
in damage accumulation between variable fatigue loading and 
constant amplitude loading to interaction effects . 

1 . Residual Stresses 

If a load amplitude is of sufficient magnitude to 
cause plastic deformation in the region of a crack tip, re- 
sidual stresses will remain in this region unless the defor- 
mation is fully reversed. These residual stresses have to 
be added to the stresses induced by the applied loads and 
as a consequence can significantly affect the fatigue damage 
accumulation . 

2 . Crack Closure Theory 

If a cracked sheet is loaded in compression the 
crack will be closed. Elber [Ref. 3] observed that crack 
closure may occur while the sheet is still loaded in tension. 
According to Elber, plastic elongation will occur in the plas- 
tic zone of the growing crack. This plastic deformation 
will remain present in the wake of the crack and will cause 
closure before complete unloading of the specimen. 

3 . A pplication 

Highly maneuverable aircraft experience a number of 
high g loadings follov/ed by a large number of smaller positive 
g loadings. Prior to completion of crack initiation these 
high g loadings produce a residual stress at points of high 
stress concentration which effectively lower the mean values 
of the subseqvient lai-gei' number of smaller positive g loads 
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and consequently reduce the actual fatigue damage . Once 
crack Initiation is complete, these high g loadings produce 
residual stress in the crack tips, which tends to retard 
the propagation of the crack and consequently reduce the 
actual fatigue damage. If an extreme negative loading 
were experienced, the residual compressive stress Induced 
would have the opposite effect and tend to cancel the favor- 
able effects from the extreme positive loading. Figures 
(la and lb) demonstrate the results Schijve [Ref. 1] obtained 
from subjecting specimens of a known Initial crack length 
to the different loading spectra. In Figure lb note the 
significant number of additional cycles (point c) needed to 
obtain the same crack length after a few high positive 
loadings. The compressive residual stresses induced by 
high loadings seems to offer a plausible reason for the 
highly conservative estimates of aircraft fatigue damage 
from a Miner's Rule approximation. It should be noted 
that the damaging ground-air-ground cycle tends to relieve 
the positive effects of extremely high maneuver load 
residual stresses and should be account for in fatigue 
damage prediction. 

In view of the effect that load sequence has upon 
fatigue prediction, it appears that not only are the loadings 
experienced by the aircraft important in damage prediction, 
but also the order in v.'hich they are imposed. A later sec- 
tion v;ill dJscuss tlie Engineering Science Data Unit (ESDU) 
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Cumulative Damage Hypothesis, which takes Into account 
sequence effects. The next section will address the prob- 
lem of estimating the actual load spectrum experienced by 
a given aircraft. 
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Ill . DATA ACQUISITION 



In transport aircraft the gust loading spectrum Is 
considered to be the most significant loading spectrum that 
contributes to fatigue damage. These loadings appear to be 
Gaussian In nature and occur symmetrically about the one-g 
level. Consequently, statistical methods such as the pov;er 
spectral density are being employed to determine the spectrum 
of loads that these aircraft have experienced. However In 
high performance military aircraft, maneuver loads are con- 
sidered the damaging loads for fatigue damage, and the gust 
loads are either disregarded or considered as a negligible 
contributor to fatigue damage. Leybold's attempt [Ref. l4] 
to describe the maneuver load spectrum from known probability 
distribution functions met with some success but v/as re- 
stricted to a particular type of aircraft. Data [Ref. 11 ] 
have shown significant differences In maneuver loads v;lthln 
aircraft of the same type as well as between aircraft of 
different types [Ref. 8]. The main reason for this differ- 
ence Is the wide range of missions performed by these air- 
craft and yet knowledge of the loads experienced by Individual 
aircraft Is necessary for accurate damage predictions. 

Stresses at fatigue critical points and their nuniber 
of occurrences are required for fatigue damage prediction. 
Stress at a point cannot be recorded directly and conse- 
quently parameters that relate to stress must bo utilised. 



The most widely used parameter for calculating stress 
Is strain. The theory of elasticity has mathematically 
developed the relationship of stress to strain. Unfortunately 
previous attempts to record strain history from strain gauges 
have been unsuccessful due to the unreliability of gauge 
readings caused by zero drift and the difficulty in processing 
and reducing strain data for a large fleet of aircraft. 

Scratch gauges have also been recommended and are capable of 
recording many hours of strain on a tiny disk but have the 
same logistics and reliability problems as strain gauges. 
Despite previous setbacks, strain recording devices are still 
considered to be the best means of obtaining stress. The 
National Aerospace Laboratory, Netherlands has developed 
a prototype system for in-flight strain recording with on- 
board data reduction [Ref. ^]. In this system the output of 
an ordinary resistance strain gauge bridge is amplified and 
digitized. This digital record is next analyzed for peak 
values using a "threshold" criterion. The values of succes- 
sive minima and maxima are written on an on-board magnetic 
tape recorder. The prototype is currently installed in an 
operational F-IQ^^G aircraft for long-term operational testing. 
Initial results are highly promising. 

Another method for detei’mlning stress, and the most widely 
used parameter for calculating stress on aircraft, is acceler- 
ation. Depending on the model of the accelerometer, discrete' 
levels of g's exceeded arcf digitaDly recorded. .Stress analy- 
sis of the fatigue critical points pi*ovldes a relationsliip 
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between acceleration and stress. This method is presently 
used by the U.S. Navy. 

A. COUNTING METHODS 

Actual load-time histories v/lll consist of a number of 
load excursions with an Irregular pattern and in an Irregular 
sequence. The analysis of load-time histories has to be 
such that the amount of damage caused by these load excursions 
is quantitatively reflected in the final results. Reference 
5 describes ten different counting procedures and their 
ability to record significant fatigue loads on aircraft. 

The methods are subdivided into three distinct categories 
depending on the parameters of the loading recorded. Basi- 
cally the counting procedures record either peaks (maximum 
or minimum values in Figure 2), level exceedances (for exam- 
ple levels one or two in Figure 2 ), Or ranges (for example 
the distance betxifeen level one and minus one). The pre- 
ferred method of counting is the range-pair-range method 
where all ranges of significance are counted. The drawback 
to this method is the extensive data reduction required. 

The restricted level-crossing count method is also con- 
sidered an effective means of recording fatigue damage. 
Accelerometers used in fatigue analysis of aircraft employ 
this method of counting. A crossing of a Devel with positive 
(or negative) slope is not made until the load has also 
crossed a second lov/er (or higher) preset level in the 
opposite direction. Intery)r'etation is hampered by 
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intermediate load cycles demonstrated in Figure 2. Two 
different load patterns are depicted but will produce equal 
counting results. However the loss of intermediate values 
is not a serious setback according to van Dijk [Ref. 5], v/ho 
concluded that the majority of the load fluctuations of 
Interest were the larger excursions that originate from the 
steady-flight level. 

B. ACCELEROMETER LIMITATIONS 

The relationship between stress and acceleration at the 
center of gravity is not as explicit as the relationship 
between stress and the strain measured at the critical 
point. The stress-acceleration relationship varies with 
weight, speed, altitude, and flight configuration. However, 
Lambert [Ref. 6] does not consider the variation critical 
for high performance aircraft since it can be assumed that 
all maneuver loads which will cause significant fatigue 
damage occur in the operational exercise of the sortie. 

Gust and maneuver loads that occur during the transit phase 
are assumed to be negligible by comparison. The stress- 
acceleration relatlonsliip is considerably simplified by a 
proper choice of an operational weight and flight parameters. 
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IV. FATIGUE MONITORING PROGRAM 



The trend to extended operational service lives of mili- 
tary aircraft emphasizes the need for accurate fatigue damage 
calculations. Since protection against fatigue does not end 
with production delivery, subsequent responsibility Includes 
further recording of service loads and environments, moni- 
toring and reevaluation In circumstances of mission changes. 
Inspection and maintenance, and structural alterations 
throughout the service life. Due to the Importance of 
fatigue strength In extending the service life of military 
aircraft, the U.S. Navy has developed an extensive fatigue 
monitoring program. V7hlch presently monitors the Individual 
structural fatigue status of 2700 aircraft. Expansion of 
the program Is planned to monitor structural fatigue life 
expended on all Navy aircraft with the exception of transport 
and rotary wing aircraft . 

The Naval Aircraft Structural Fatigue Life Program enable 
structural Integrity decisions to be made realistically In 
terms of fatigue experienced by the structure rather than 
by continuing the previously nebulous procedure of using 
the number of years in service or the accumulation of flight 
hours . as an expended life criterion . The Naval Air Develop- 
ment Center Is tasked with monitoring the structural fatigue 
life expended on eacli Naval aircraft. As a result a fleet- 
wide counting .accelerometer program, combined with pilot's 



reports relating flight usage data, laboratory fatigue test 
data, and cumulative fatigue damage theory, was developed 
to provide a measure of the actual service life expended for an 
Individual aircraft throughout Its service life. 

A. UTILIZATION 

Usage of the fatigue data accumulated has provided the 
Navy with the source of information necessary to effectively 
manage a large fleet of aircraft. The data are reported 
quarterly [Ref. 6] and used for: 

Early recognition of changes in service usage trends. 

Identification and monitoring of aircraft accumulating 
unusually severe or excessive load factor occurrences and 
aircraft v;ith unusually high fatigue damage accumulation 
rates . 

Determination of aircraft time of retirement or 
rotation into or out of severe service usage . 

Extension of service life based on changes in 
mission, flight restrictions, replacement or modification 
of major structural components, review of service history, 
and performance of full scale fatigue tests. 

Provision of information effecting both organisational 
and in-service maintenance problems. 

Utilization of data in the development/modification 
of naval aircraft design and test specif ications/requlrements . 
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B. DAMAGE PREDICTION 



As was previously mentioned, the main function of the 
structural fatigue program is to monitor the damage accumu- 
lation of individual aircraft. Typical of the method used 
to arrive at the damage calculation is the method employed 
to calculate the damage on A- 7 B aircraft. The following 
discussion will specifically apply to A- 7 B aircraft but is 
considered representative of the other types of aircraft 
monitored under the program. 

1 , Full Scale Testing 

The first step in developing a damage calculation 
on the A- 7 B was to determine fatigue critical points. These 
points were analyzed and S~N data appropriate to these points 
were developed. A block program based on estimated loading 
per 100 hours of life time was used to fatigue test the 
aircraft (Table I). 



TABLE I 

A -7 Block Loading 



I ^ 

Maximum Number of load applications 

Load (g’s) per load block 





Block 1 -^ 


Block 5-10 


2.^5 


1700 


1700 


3.15 


950 


950 


3.85 


650 


650 


^.55 


^50 


^150 


5.25 


250 


250 


5.95 


] 36 


136 


6.65 






7.35 


15 


15 


8.05 






8.75 


]. 


2 






Based on the proposed block loading pi’ogram an 
estimated ' fatigue life was predicted at each critical point 
based on Miner’s cumulative damage theory. The full-scale 
fatigue test was then conducted (using blocks 1-10, which 
represent 1000 hours of flight time) until a crack initiated 
at one of the critical points. The stress concentration 
factor of the S-N data was then adjusted to predict a fatigue 
life of 100 percent for the number of cycles that occurred 
on the full-scale fatigue test. The loads obtained from 
data reduction of accelerometer readings are then applied 
to the S-N data and a cumulative damage prediction is calcu- 
lated. Due to the uncertainty of the actual flight loading 
history a safety factor of two is applied to the dam.age 
prediction . 

2 . Load Data Reduction 

The accelerometers used on the aircraft record the 
total exceedances at only four preset levels. The custodian 
squadrons report quarterly to NADC the accelerometer readings 
and pertinent flight data. Based on the four exceedance 
levels recorded, seven other acceleration levels are computed 
for the loading spectrum to be used for fatigue prediction 
in the case of A-7 aircraft. The computed exceedances are 
determined by use of the following equation [Ref. 8'J: 

Cumulative Counts = Ae^'^^ (3) 
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Table II provides the values of the constants to be used 
for the seven levels. 



TABLE II 

A-7 Load Level Approximation 



Cumulative Counts 


A 




B 


C 


5.5 


g' s 


1.07 


ln(5-g 


counts) tin (6-g counts) 




6.5 


g’ s 


1.00 


ln(6-g 


counts)+ln(7-g counts) 


2 


7.5 


g’s 


1.00 


ln(7-g 


counts ) +ln( 8-g counts) 


2 


8.5 


g' S 


1.00 


2 In ( 8-g 


counts)-ln(7 .5-g counts) 


1 


^.5 


g' S 


.82 


31n(5-g 


counts ) “ln( 6-g counts) 


2 


^.0 


g' s 


.97 


21n(4 .5- 


-g counts )-ln(5“g counts) 


1 


3.5 


g' S 


.95 


21n ( 4-g 


counts )-ln( . 5-g counts) 


1 



The basis for these equations is obtained from reduced 
oscillograph data. An oscillograph recorder system is 
Installed on a number of fleet aircraft which records time 
histories of airspeed, altitude, and normal acceleration. 

The recorder is energized by a switch on the landing gear 
retraction mechanism and consequently records data only v;hen 
the landing gear is in the up position (the same is also 
true of accelerometer data so no landing information is pro- 
vided) . Periodically these oscillograph data are compiled 
on aircraft [Ref. 8] and used to develop relationships betv/een 
the various levels. The data .are also used to redefine the 
stress-acceleration relationship since altitude and airspeed 
data are correlated v/lth the loads . Unfortunately a semi- 
automatic data reduction system is employed which outputs 
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data on all A-7’s, and consequently, information on indivi- 
dual aircraft is not available for a comparative analysis of 
trends betv/een individual aircraft. The 11 levels of accel- 
eration now available are transformed into stress levels 
and the S-N data ,are referred to for damage prediction. The 
entire process has been adapted to a computer system which 
enables rapid retrieval of pertinent flight and fatigue data 
of any aircraft currently under this program. 

C . IMPROVEMENTS 

The fatigue monitoring program conducted by NADC has 
unquestionably made significant advancements towards an 
accurate prediction of fatigue damage. There are possibili- 
ties for improvement based upon current knowledge of the 
fatigue phenomenon. 

It has been determined that maneuver loads occur randomly 
and that block testing will not produce the same fatigue 
damage as random load testing. Schljve [Ref. 1] mentions 
that, in general, the life in "equivalent" program tests 
is longer than in random tests. He alludes to one series 
of programs, v.’here fatigue lives v^ere about six times longer 
than in random load tests . 

Through block loading fatigue testing on an A-7 aircraft. 
Miner's Rule has been altered by shafting the S-N data to 
account, somewhat, for the variation of load amplitude 
throughout the flight profile. Where Miner's Rule may be 
too conservative, block loading on the other hand may be 
too optimistic. 



The sequence effects of random loading cannot be 
ignored if accurate predictions are to be obtained. The 
remainder of this thesis considers the feasibility of 
applying sequential effects for fatigue analysis based 
on four-level accelerometer readings and the suitability 
of applying this method to a computer. The next section 
will address the problem of developing a reasonable flight 
profile suitable for application to sequence effects and 
based on four levels of acceleration provided. 



V. DEVELOPMENT OF A PLIGHT PROFILE 



In order to apply sequential effects a reasonable flight 
profile must be produced that satisfies the following 
requirements : 

1) The order In which the loads are encountered must 
be randomly distributed. 

2) The loading spectrum must be generated from the 
four levels provided by the accelerometer readings. 

3) Recognizing the Individuality of aircraft flight 
profiles, two equivalent sets of accelerometer 
readings must not necessarily generate the same 
loading spectrum. 

Reference [8] and other similar reports from NADC con- 
taining oscillograph Information on aircraft were analyzed 
In developing a flight profile suitable for high performance 
military aircraft . 

A. OSCILLOGRAPH/ACCELEROMETER DATA 

Schljve [Ref. 1] In Figure 3 differentiated between the 
gust spectrum and the maneuver spectrum experienced by air- 
craft. The gust load spectrum consists of an equal number 
of positive and negative small amplitude peaks occurlng about 
the steady flight condition. The maneuver load spectrum 
consists of the remaining loads an aircraft experiences and 
Is highly asymmetrical about the steady flight condition. 
Lambert [Ref. 6] made a similar analogy as shown In Figure ^4 



?6 



and attributed the negative , loads In the maneuver spectrum 
as being due to "overswlng" resulting from recovery from a 
hlgh-g maneuver. Prom an extensive study of Swiss VENOM 
aircraft Branger [Ref. 9] concluded that the number of 
maximum peaks above the steady flight condition was equal 
to the number of minimum peaks belov7 the steady flight 
condition but that the distribution Is highly asymmetric. 
Oscillograph data [Ref. 8] support this phenomenon. A 
lower threshold of two g's in the positive loading regime 
and an upper threshold level of .5 g's in the negative 
flight regime produced ratios of positive tq negative g's 
as high as eight to one. The reason for this is demonstrated 
in Figure 3 where the number of positive loadings above 
threshold level A clearly exceed the number of negative 
loadings below threshold level B. Schijve [Ref. 1] indicates 
that all significant fatigue damaging loads originate from 
around the one-g level . Based on this evidence all 
accelerometer readings are assumed to be maximum peaks and 
llkevjlse all oscillograph peaks in the positive loading 
regime are also assumed to be maximum peaks since the 
criterion for recording a peak is essentially the same 
as the accelerometer criterion. All loadings are assum.ed 
to originate from the one-g level or less (minimum peaks) 
depending on the number of negative regime peaks the 
oscillograph data provides. 
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B. ANALYSIS OF DATA 



Different assumptions are made relating the peaks to a 
common trend and the data are analyzed to verify the 
validity of each assumption. 

1 . Cumulative Loads Between Types 

Since all the aircraft considered were high perform- 
ance aircraft, having missions that required high maneuvering 
loads, it v/as assumed that all aircraft experienced essen- 
tially the same number of loads based on 1000 hours of flight 
time. Figure 5 indicates the logarithm of the number of 
loads that exceeded a particular g level for A-AE, A-7A, 

F-i<B, A-6A and F-8E aircraft per 1000 hours of flight time. 
The data would not support the assumption that the maneuver 
loads of each type of aircraft were essentially the same and 
thus the mission requirements of each type of aircraft are 
not comparable enough to produce an equal number of similar 
loads. Hov/ever, due to the similarity of the curves the 
relative ratios of load levels between different types of 
aircraft may be proportional. 

2 . Ratio of Loads Between Types 

Although all the aircraft considered had m.issions 
demanding high maneuver loads, some aircraft had missions 
that demanded more loads per flight than those of another 
aircraft type. In spite of this it was assumed the loading 
spectra were similar and the ratios of loading levels 
botv;een the types of aircraft are essentially equal. All 
loading levels \:ere ratloed to the lowest posiltlve load 
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level recorded. Figure 6 shows the ratio of five load 
levels for A-^E, A-6A, A--7B, F-^B, and F-8E to their respec- 
tive 2. 0-2. 5 g load counts. Two of the aircraft’s ratios 
deviated considerably from the other three aircraft’s ratios. 
Intermediate load levels were then investigated to avoid the 
possible effect of extreme gust loads at the lower level. 
Figure 7 shows the ratio of four load levels to the 4-4.5 g 
load level. The ratios again Indicated a considerable degree 
of scatter. All loads that exceeded an upper level were 
ratloed to the load counts at lower levels. Figure 8 shows 
the ratio of four load levels to the number of loads that 
exceeded the four g load level. These ratios indicated even 
more scatter than those obtained in Figure 7- 

Since the load latios investigated did not indicate 
any clear similarity , it was concluded that the flight 
profiles between types of aircraft were distinctly different. 

3 . Extreme Values 

Extreme values were investigated to determine whether 
there was a relationship between the load exceeded by extreme 
positive and negative loads. It was assumed that entry into 
these extreme loads vjas not intentional and probably flights 
that encountered these loads encountered extreme vaD.ues in 
both directions. Extreme values were defined as loads that 
accounted for one percent of the total number of loads 
occurring in either the positive or negative regime. Table 
III lists the number of positive and negative levels exceeded 
one percent of the time. 
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TABLE III 



Extreme Value Occurences and Level of Exceedance 



A/C Type 


Level 


No. Exceed. 


Level 


No . Ex 


A-4E 


6.0 


282.6 


-.41 


17.2 


A-6A 


5.5 


88.9 


-.21 


4.9 


A-7A 


5.5 


177.5 


-.1 


20.3 


P-4B 


5.5 


1 — 1 
C\J 


-.41 


8.1 


F-8E 


6.0 


117.8 


-.41 


35.4 



Analysis of Table III indicates no definite relation- 
ship between the level exceeded in the positive regime and 
the level exceeded in the negative regime. A sample corre- 
lation coefficient of .642 for levels of extreme-value 
exceedance was not considered conclusive enough to support 
the assumption that a relationship existed between the 
positive and negative levels. Analysis of extreme values 
on a flight-by-fllght basis is needed to confirm the 
assumption . 

Reference 10 contains oscillograph data on the Navy’s 
"Blue Angels" Plight Demonstration Team compiled on a flight- 
by-flight basis. These data were also analyzed for a possible 
trend. By analyzing the data on an individual flight basis 
it could be determined whether or not extreme high loads 
tend to induce corresponding extreme low loads as a result 
of over-correcting recovery from the high load maneuvers. 

TabDe IV lists the highest and lowest g level recoi'ded for 
12 flights analyzed at random. 
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TABLE IV 



Blue Angel's Flight- 

Flight 

1 

2 

3 

5 

6 

7 

8 

9 

10 

11 

12 



y_Flight Extreme Values 



Max . g 


Min . g 


7.^ 


- 1.56 


6.18 


- 1.05 


5.83 


- 1.^1 


i|.6 


- 1.56 


5.62 


- 2.87 


6.il8 


- .11 


6.i|8 


- .05 


6.15 


.0^ 


7.38 


- 1.27 


6.63 


.08 


5.22 


- 1.15 


e.i\9 


- 1.91 



Detailed analysis of Reference 10 Indicated that extreme low 
loads neither preceded nor followed an extremely high load 
and as Table IV indicates there is absolutely no foundation 
for assuming a high probability of both high and lov-; values 
occurred on the same flight. Of particular significance in 
Table IV is the fact that the lowest negative reading (-2.87) 
occurred during a flight v/hich had the third lovjest positive 
reading. It must be concluded that the negative readings 
occur randomly and independently of any extreme high loading 
conditions; consequently, they could be considered as part 
of a "basic" Heading spectrum unJque to the flying qualities 
of the particular aircraft consido’.’cd . Tlie oxti’ome How g 
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loads are probably due to high gusts loads and therefore 
the gust loads should be considered Integrated Into the 
'•basic*' loading spectrum. The dynamic "overswlng" Lambert 
[Ref. 6] discussed must be above the .5 g level threshold. 
From the data investigated so far it appears that each air- 
craft has a flight loading profile distinct from another 
type of aircraft and this profile Includes all the negative 
g loadings and excludes all high g loads. If types of 
aircraft do Indeed have a loading spectrum "fingerprint", 
it can only be determined by comparing oscillograph data of 
individual aircraft of the same type . 

4 . Individual Aircraft by Type 

In order to substantiate whether aircraft of the 
same type have identical flight profiles below some extreme 
positive load level, individual aircraft dataware needed. 

A study conducted for NADC by Technology Incorporated [Ref. 

11] contains Individual aircraft oscillograph data on four 
A-6A's and six A-7B's collected over a one-year period and 
representing data on 306.9 flight hours for the A-6 A's 
and 501.9 flight hours for the A-7B's. Although the number 
of aircraft of one type analyzed was small it was considered 
significant enough to determine if aircraft of the same 
type do experience essentially the same loading profile 
below some high load level. 

a. Validity of Data 

The data were .first compared v/.1th the osclllograi h 
data of Ref. 8 to assure comparable results betv.'eon the two 
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reports. The A-6A data agreed quite closely and the A-7B 
data also agreed favorably In the lower loading regime but 
diverged significantly above the flve-g level. This was 
not considered Important In this Investigation, since the 
region of Interest for developing a loading profile Is 
contained below the flve-g level. Since the ratios of loads 
were considered to be the key to developing a flight profile, 
the data were considered to agree favorably with other 
oscillograph data of the same type aircraft . 

- b . Comparisons of Each Type 

The number of exceedances for Individual aircraft 
of both types at each level were compared. The exceedance 
distribution for Individual aircraft of each type Indicated 
a high degree of similarity which would be necessary to 
develop a flight profile for each type. Figure 9 shoves the 
number of loads exceeding an indicated level per 10,00 hours 
of flight tlm.e for each of the A-7B aircraft. The negative 
loads also exhibited a comparable similarity. Analogous 
results were also obtained for the four A-6A aircraft. 

The number of loads occurring In each range vjere 
ratloed to an arbitrary level that had a range of one g. 

The one-g range level was chosen to agree v;lth the ranges 
bet\-jeen recording levels of accelerometers installed In the 
aircraft. The number of loads that occurred within a prescribed 
range v/ere determined by subtracting the number of exceedances 
at the higher level of the range from the number of exceedances 
at the lov;er level of the I’ange . The~one g range level wa.5; 
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increased and comparisons made until an unacceptable dis- 
agreement resulted between individual aircraft . Acceptable 
agreement for the A-6A occurred only for readings up to the 

3- ^ g range. However acceptable agreement for the A-7B 
occurred up to the 4-5 g range. Figure 10 shows the ratio 
of load occurrences to the number of 4-5 g occurrences for 
the A-7B aircraft. The load values vjere represented by the 
mean value of each range. The fact that one aircraft exper- 
ienced an increase in the ratio from the 2. 5-2. 9 range to 
the 2. 9-3. 3 range was not considered significant since that 
aircraft had recorded only 19.7 hours of flight time with 
the oscillograph installed. It was considered significant 
that four of the six aircraft had almost identical values. 

c. Plight Profile 

Based on the above results it v;as decided to 
develop a flight profile for A-7B aircraft that would depend 
solely on the accelerometer counts in the 4-5 range. Table V 
Indicates the load range and the ratio of that range to the 

4- 5 g’s each aircraft experienced. The ratio of loads in 
each range were assumed to be normally distributed about the 
mean ratio. Based on these data a mean value and a sample 
standard deviation were calculated. It was decided that a 
range of 1.65 times the sample standard deviation from each 
ratio mean would encompass an adequate range of loads that 
the aircraft would experience in each load level. Multiplying 
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TABLE V 



A- 7 B 


Individual 


Ratios 


of 


Positive Loads 


to 4-5 g Loads 


Load 


19.7 


50.2 




140.5 


101.3 


113.1 


76.9 


Range 


hrs . 


hrs . 




hrs . 


hrs . 


hrs . 


hrs . 


2 . 0 - 2. 5 


3.74 


3.81 




3.29 


5.24 


5.14 


2.85 


2 . 5 - 2. 9 


1.87 


1.67 




1.86 


2.67 


1.76 


1.73 


2. 9 - 3 . 3 


1.99 


1.18 




1.46 


2.51 


1.27 


1.50 


3. 3-3.6 


1.14 


.658 




.823 


1.49 


.714 


.855 


3. 6-3. 9 


.73 


.504 




.621 


.770 


.559 


.605 


3. 9-4. 2 


.444 


.372 




.438 


.414 


.400 


.439 


4. 2-4. 5 


' .298 


.333 




.339 


.304 


.347 


.325 


4. 5-4. 8 


.275 


.268 




.252 


.275 


.257 


.251 


4. 8-5.0 


.152 


.149 




.143 


.181 


.143 


.135 


each ratio by the 


number 


of 


4-5 g 


readings 


would produce 


the 



actual number of loads the aircraft experienced for that 
flight period. 3.3 times the sample standard deviation pro- 
vides a probability of .95 that the aircraft experienced a 
load level in that range based on a sample size of six. 

Table VI gives the range of ratios for each load level that 
the aircraft is allowed to experience in the positive regime 
for each 4-5 S reading. 

Table VII gives a similar list of ratios as Table 
V for the negative regime. The excessive number of negative 
readings aircraft (Serial No. 154370) experienced were not 
considered in the calculation of mean and sample standard 
deviation of the negative regime ratios. This aircraft had 
only 19.7 hours of recorded load time and its ratios v/ere 
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TABLE VI 



A- 7 B Positive Ratio Range For Each Load Level 



Range Mean 


Low 


Mean 


High 


2.25 


2.399 


.01 


5.62 


2.7 


1.306 


1.92 


2.53 


3.1 


.815 


1.65 


2.48 


3.^5 


.429 


. 9^7 


1.47 


3.75 


.il 65 


.6315 


00 




.370 


.417 


.464 


^.35 


.292 


.324 


.356 


k.65 




on 

C\J 


.281 


^.9 


.122 


.1505 


.179 



A- 7 B Individual 


TABLE VII 
Ratios of Negative 


Load Range 


to 4-5 g 


Load 


Load 


19.7 


50.2 


140.5 


101.3 


113.1 


76.9 


Range 


hrs . 


hrs . 


hrs . 


hrs . 


hrs . 


hrs . 


. 6 -. 4 


18.97 


9.02 


4.807 


9.24 


6.74 


5.97 


. 4-. 2 


4.66 


2.96 


1.506 


1.88 


1.412 


.606 


. 2-.0 


1.22 


.601 


.376 


.408 


.329 


.080 


0 -(-. 2 ) 


.274 


.175 


.137 


.110 


.071 


.020 



not considered representative of the typical A- 7 B. Neglecting 
the excessive number of negative loads on this aircraft would 
not be more damaging in a cumulative damage calculation since 
• sample calculations with representative S--N data indicate a 
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loading with a larger range but a lower mean was less 
damaging than a loading with the same maximum peak that 
originated at one g and had a smaller range. Table VIII 
Indicates the ranges of ratios in the negative regime that 
have a probability of .95 of occurring for each 4-5 g 
reading based on a sample size of 6. 

TABLE VIII 

A-7B Negative Ratio Range For Each Load Level 



Range Value 


Low 


Mean 


High 


.538 


3.97 


7.155 


10.34 


.346 


.260 


1.673 


3.09 


.148 


.050 


.359 


.667 


.054 


.004 


.103 


.200 



Figures 11 and 12 display the range limits and 
mean for each load level for both the positive and negative 
regime, respectively . The mean value of each positive load 
level Is considered to be an adequate representation of the 
loads in that range since the greater number of loads that 
occur below the mean value are offset in the damage calcula- 
tion by the fewer but more damaging high values In each 
range . 

5 . Sequence of Loads 

Based, on the previous constraints a computer program 
was developed to arbitrarily assign a I'atlo for each load 
level to the 4-5 g level reading. The value of each ratio 
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was determined by means of a normally distributed random 
number generator. The number of actual loads In each load 
level v;ere then determined by multiplying each ratio by 
the number of 4-5 g readings that were Inputted. 

Three other levels of loads from accelerometer 
readings were also inputted that represented loads that 
occurred in the 5-6, 6-7, and 7-8 load ranges respectively. 
These load ranges were arbitrarily assigned mean values of 
5.75, 6.6, and 7.8 which v;ere considered high enough to 
represent a conservative estimate for a damage calculation 
of a large number of loads occurring in this range since 
75 percent of loads recorded occurred below these values 
in each range [Ref. 11]. No attempt was made to distribute 
these high level loads throughout the entire range of their 
possible values due to the limited number of occurrences 
of these loads in Ref. 11. However if distribution functions 
for the extreme value loads were known it would be appropriate 
to subdivide the loads occurring in each range into a smaller 
number of load ranges. Distributions similar to the one 
used in Table II could be used but in this particular case 
the distribution could not be substantiated from the data 
in Ref. 11. 

Based on the four load level accelerometer readings, 
all loads in the positive regime were assigned to be maximum 
peaks. Associated vjlth the maximum peaks v;as an equal number 
of minimum peaks v;hlch vjere assigned values from the negative 
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regime and the remaining number of required minima were 
assigned a value of one g. 

A uniformly distributed random number generator 
placed each group of peaks in a random order. Starting 
with a minimum peak, a load value V'/as picked from each group 
until all the loads in each group had been chosen. A final 
half cycle was then assigned to finish at a one-g level. 
Figures 13-16 indicate four flight profiles generated for 
accelerometer readings of three, five, three, and four 
representing loads occurring in the 4-5, 5-6, 6-7, and 7-8 
g ranges respectively. 

It is highly unlikely that the order of loads gener- 
ated from the accelerometer readings would represent the 
actual sequence the aircraft experienced during that period 
and thus a damage calculation based on this loading profile 
would be inaccurate. However, as the aircraft experiences 
an increased number of loadings over its fatigue life it is 
highly probable that the aircraft experiences a majority 
of the sequential loadings generated or ones quite similar 
to them since the loads generated are based on actual data 
obtained from aircraft of the same type whose order of 
loading is known to occur randomly. The flight profiles 
generated will be of increasing significance during the 
latter stages of the aircraft ’ s ' life where the accuracy 
of a damage prediction calculation becomes of increasingly 
higher importance . 
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6 . Suitability of Flight Profile 



The requirements stated previously for generation 
of a profile have been met. The loads are randomly distri- 
buted and are generated- solely from accelerometer readings. 

Two equivalent accelerometer readings do not necessarily 
generate the same loading profile. Due to the small number 
of aircraft sampled it cannot be assumed that the above pro- 
file is an accurate representation of the majority of A-7B's 
profiles. A larger sampling of aircraft would be necessary 
to further refine a representative profile. Currently the 
lowest accelerometer level being r.ecorded on A-7B’s is the 
flve-g level. In order to adapt the flight profile developed, 
the lowest level recorded would have to be lowered to the 
four-g level. The same is also true of the A-6A where the 
lowest level presently recorded is the four-g level and a 
three-g level is considered necessary to obtain an indicator 
of the A-6A profile. However the above procedure demonstrates 
the feasibility of generating a flight profile from acceler- 
ometer readings since each type of aircraft apparently does 
have a "fingerprint" load spectrum below some high load 
level. The flight profile developed above will be considered 
representative enough of high performance fleet aircraft to 
proceed vjith accounting for the sequential effects in a 
fatigue damage calculation. The next section deals with 
the feasibility of applying a fatigue damage prediction that 
accounts for sequential effects and is suitable for computer 
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VI. ESDU CUMULATIVE DAMAGE HYPOTHESIS 



A cumulative damage hypothesis is needed that accounts 
for the order in which loads are applied to an aircraft that 
is suitable for computer analysis. 

The ESDU Cumulative Damage Hypothesis [Ref. 2] meets 
these requirements. Although originally intended for block 
loadings the method has been modified for analysis of loading 
on a cycle by cycle basis. The ESDU method requires no 
additional data over that required for use in Miner’s Rule. 
Testing [Ref. 2] of lab specimens under block program loading 
has produced results that demonstrated a far better predic- 
tion of fatigue failure by the ESDU method than by Miner's 
Rule. A number of similar specimens were tested until failure 
which on the average took 30.9 programs to occur. i The Miner's 
Rule estim.ate for failure was 11.7 programs; whereas the ESDU 
estimate was 28.7 programs, an improvement of over 2^0 
percent . 

The basic hypothesis of the ESDU method is centered in 
the concept that, once the stress applied to a component ex- 
ceeds the proportional limit, residual stresses build up 
vjhich effectively alter the actual mean stress of subsequent 
stress cycles. The following basic assumptions are made 
when applying the ESDU method: 

1) The material is perfectly elastic belov/ the 
yield stress and perfectly plastic above it. The yield 
stress is taken as the .2 percent offset stress. 



2) The stress-strain characteristics of the material 
are the same in tension and compression. 

3) The stress concentration is constant for all 
nominal stresses. 

\ 

Yielding is localized to stress concentrations 
and therefore does not significantly affect stresses in the 
rest of the member. 

5) A constant amplitude S-N curve for the material 
is available. 

A member that is subjected to a nominal stress whose 
peak value is SMAX has an elastic stress at the stress 
concentration of SMAXEL equal to: 

SMAXEL = KT • (SMAX) (3) 

If the value of SMAXEL exceeds PP, the .2 percent offset 
yield stress, the mean value of the subsequent cycles will 
be lowered by an amount equal to H where: 

H = FP - SMAXEL (4) 

The same is also true for an excessive compressive stress, 
which due to the build up of residual tensile, stress causes 
an effective raising of the mean stresses of the subsequent 
cycles. Figures 17--19 indicate the resulting effect that 
the different types of loadings have on subsequent cycles. 

In Figure 17 neither the maximum nor the minimum peak exceed 



PP so no alteration of the mean stress in subsequent cycles 
is required. Figure l8 indicates the effective raising of 
the mean stress of subsequent cycles due to a compressive 
stress that exceeded -FP . Figure 19 shows the effect of 
compressive residual stresses due to a tensile loading 
whose maximum value exceeded +PP. 

Application of the ESDU method is quite simple and is 
easily adapted to a cycle-by-cycle analysis. A possible 
objection to the validity may be that the effects of stress 
relaxation are neglected, which causes the residual stress 
to fade avjay if given enough' time. However Schljve [Ref. 1] 
states that this effect has not been observed in aluminum 
alloys, the materials most aircraft wing components are 
made from. 

The next section discusses the computer program used to 
develop a flight profile from accelerometer readings and to 
apply the ESDU Cumulative Damage Hypothesis for a damage 
calculation of this profile. 



VII. DISCUSSION OF COMPUTER PROGRAM 



Program ESDU was developed primarily to demonstrate the 
feasibility of adapting the ESDU Cumulative Damage Hypothesis 
to a computer program for use in calculating the fatigue 
damage of high performance aircraft with accelerometer 
readings as data input. In order to provide a concrete 
example the program v/as specifically designed to provide 
damage calculations on A-JB aircraft that had accelerometer 
recording levels at the four, five, six, and seven g level. 

A. FEATURES 

There is virtually no limitation to the number of readings 
at each level (ten digit maximum); however, only a limited 
number of readings will be used in the calculation if the 
number of four-g level readings exceeds 325. The program 
requires 152 K bytes of storage and less than three seconds 
of computer time for each aircraft calculation. Data input 
consists of the number of aircraft to be analyzed, the 
serial number and four accelerometer readings for each 
aircraft. 

B. PROCEDURES FOR CALCULATION 

The program is essentially subroutine structured with 
the main program linking the ten subroutines required for 
the calculation. The main program accepts the input data, 
assures compatible storage requirements, requests maximum 



and minimum peaks, requests a random ordering of the peaks, 
and finally calls for a calculation of damage for the flight 
sequence using the ESDU method. A brief description of 
each subroutine follows . 

1 . Subroutine Fix 

Fix assures the storage requirements will not exceed 
152 K bytes. If more than 325 four-g level readings are 
Inputted for one aircraft, it is hypothetically possible to 
generate more than the 10,000 peak storage limitation. If 
more than 325 four-g level loads are inputted the four 
recorded levels are ratloed so that 325 four-g level loads 
will be provided to the main program. Fix also provides a 
correction factor to compensate for the reduced number of 
peaks in the final damage calculation. 

2 . Subroutine Upspec 

Upspec provides the main program with the number of 
maximum peaks and their values to be used in the damage 
calculation. A ratio of each load level below four g’s to 
the number of 4-5 g loads is determined by a normally 
distributed random number generator provided by subroutine 
Gauss. If the random number provided exceeds the .95 
percentile range for each load level a new number is 
requested. After the number of loads in each level are 
determined, subroutine Hiload assigns the mean g load value 
to each peak in each load range . 



3 . Subroutine Dnspec 



Dnspec performs essentially the same function as 
Upspec, except it provides the main program with minimum 
peaks. Appropriate ratios are assigned to four levels of 
minimum peaks . After determination of the number of peaks 
required in each load level the remaining required minimum 
peaks are assigned to a fifth level, which represents the 
one g minima, to equalize the number of maxima and minima. 
Subroutine Dnload assigns the load value to each of the 
five ranges . 

^ . Subroutine Mixer 

Mixer provides the main program with the maxima and 
minima arranged in a random order so as to represent a flight 
profile. Mixer requests subroutine Shufle to random^ly 
arrange both the maxima and minima by means of a uniformly 
distributed random number generator provided by subroutine 
Random. Mixer then alternately picks a value from both 
groups and stores them in an array until all peaks have been 
used . 

5 . Subroutine ESDU 

ESDU calculates the damage Incurred by the aircraft 
on a cycle-by-cycle basis using the ESDU Cumulative Damage 
Hypothesis. ESDU requires data about the particular aircraft 
being analyzed. A stress concentration factor, S-N data, 
ultimate tensile stress, .2 percent offset yield stress, 
and conversion from g to stress are needed for the aircraft 
at the fatigue critical point. Polynomial approximations 
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were used for the S-N data but more precise descriptions 

of the S-N data would be desirable. Linear interpolation 

of tabulated S-N data would be an acceptable method. S-N 

data for only one mean value are all that is necessary 

.since theT>rogram adjusts' the S-N data to accommodate 

different means by the method described in Ref. 2. In the 

✓ 

next section sample data runs are discussed and the 
relevancy of their- results to actual aircraft predictions. 



VIII. DATA ANALYSIS 



Actual accelerometer readings reported from fleet air- 
craft during heavy maneuver periods were fed into the program 
to determine if reasonable damage predictions would result. 
Although it was realized that the damage predictions have 
no realistic value except for A-7B’s due to the flight pro- 
file variations between types of aircraft, subroutine ESDU 
was altered to accommodate the pertinent load data for A-6A 
and F-^B aircraft as well as A-7B aircraft. A fourth "hypo- 
thetical" aircraft was devised in order to cover a range of 
stress concentration factors and g-to-stress conversions 
commonly found in high performance aircraft. The stress 
concentration factors and g-to-stress conversions represent 
the actual values for the fatigue critical points of the 
aircraft types tested. The g-to-stress conversion is deter- 
mined by the design limit load of the aircraft, which is 
the amount of nominal stress the fatigue critical component 
experiences at the 6.5 or 7.0 g load depending upon the 
aircraft in question. The design limit loads varied from 
^0,000 p.s.i. on the A-6A to 2^,600 p.s.l. on the A-7B. 

S-N data was obtained by interpolation of Goodman Diagrams . 
contained in Ref. 12. The material used at the fatigue 
critical points for all aircraft was aluminum plate 7075-T6. 
Table IX provides the pertinent Information on each aircraft 



used in the test runs. 



TABLE IX 





Aircraft Tested 






Aircraft 


P- 4 B 


A- 7 B 


A- 6 A 


Hypo . 


KT 


3.3 


3.5 


3.0 


2.5 


stress (KSI) 
per g 


4.723 


3.5143 


6.154 


4.57 


Location 


outer 


outer 


inner 




critical pt . 


wing 


wing 


wing 


- 




panel 


panel 


panel 




Plight hours 


25.8 


31.5 


33.^ 


31.5 



Five test runs were performed on each aircraft to deter- 
mine if the prediction variation was significant with differ- 
ent flight profiles. Data for the 4-5 g level of the A- 7 B 
were developed from the equations provided in Table III . 

Each identical flight profile was subjected to a Miner’s 
Rule approxlm.ation for comparison with the ESDU method. 

Tables X-XIII show the results of these runs for each aircraft. 



TABLE X 



P- 


• 4 B Test Runs 


with 


25.8 Hours 


of Combat 


Plight Time 




Acceleration 


data 


(g’ s ) 


Damage 


prediction 


Four 


Five 


Six 


Seven 


ESDU 


Miner ' s 


30 


11 


2 


0 


.0005824 


.0156074 


30 


11 


2 


0 


.0006568 


.0153307 


30 


11 


2 


0 


.0008384 


.0157427 


30 


11 


2 


0 


.0009822 


.0157232 


30 


11 


2 


0 


.0010484 


.0155861 
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TABLB; XI 



A- 6 A Test Runs with 33 *^ Hours of Combat Plight Time 



Four 


Accelerometer 
Five Six 


.data 

Seven 


Damage 

ESDU 


prediction 
Miner ' s 


26 


5 


0 


0 


.0010131 


. 01 ^ 47^39 


26 


5 


0 


0 


.0012068 


. 01^47753 


26 


5 


0 


0 


.0011503 


.0165619 


26 


5 


0 


0 


. 00120^12 


.0158776 


26 


5 


0 


0 


. 001 ^ 105 ^ 


.0166387 









TABLE XII 






A- 7 B 


Test 


Runs with 


31.5 Hours 


of Combat 


Flight Time 




Accelerometer 


data 


Damage prediction 


Four 


Five 


Six 


Seven 


ESDU 


Miner ' s 


194 


120 


57 


12 


.00182^18 


.0118987 


194 


120 


57 


12 


.001703^ 


.0116276 


194 


120 


‘57 


12 


.0017585 


.0117723 


19^ 


120 


57 


12 


.0017833 


.0119610 


194 


120 


57 


12 


.0016321 


.0116316 
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TABLE XIII 



Hypothetical Aircraft with A- 7 B Flight Data 





Acceleration 


data 


in g' s 


Damage 


prediction 


Four 


Five 


Six 


Seven 


ESDU 


Miner' 


19 ^ 


120 


57 


12 


.0016852 


.0058^165 


19 ^ 


120 


57 


12 


.001^576 


.0059167 


19 ^ 


120 


57 


12 


.001^1858 


.0061357 


19 ^ 


120 


57 


12 


. 00136^11 


.0059771 


19 ^ 


120 


57 


12 


. 0012^137 


.0058356 



In order to easily Interpret the degree of variation 
between the five flight profiles, the mean value of the 
damage calculation and the percent of deviation from the 



mean of 


the largest 


excursion 


for 


each aircraft was 


calcu- 


lated . 


The results 


are tabulated 


in Table 


XIV. 








TABLE 


XIV 










Results of 


Test 


Runs 








Aircraft 


F-i|B 




A~ 7 B 


A- 6 A 


Hypo. 


ESDU mean 


.0008216 


• 


0017^0^1 


.001196 


. 001 ^ 4^173 


Percent 


excursion 


29.11 




. 8i| 7 


17.5 


16 .'Hi 


Miner ' s 


mean 


.015598 


• 


011778 


.015719 


. 0059 'I 23 


Percent 


excursion 


: 1.71 




1.55 


6 . 2 


3.25 



Some definite conclusion^ were made from the test runs. 
The ESDU method makes a considerable difference in a damage 
prediction as compared to Miner’s rule. The variation from 
the mean in the damage prediction was noticeably larger in 
the ESDU method than from Miner's calculation, although they 
both vrere based on the same flight profile. Apparently the 
variation in the number of loads in the profile is not quite 
as significant in the damage calculation as is the order of 
the loading since the ESDU method, which accounts for the 
order of the loads, had a much larger deviation. 

A. LIFE ESTIMATE 

It was decided to calculate the life expectancy of an 
A-7B based on the number of occurrences of loads per 1000 
hours from the oscillograph data of Ref. 8. Due to the 
large number of loads involved the data were subdivided 
into blocks of 26 representing 38.46 flight hours each. The 
loads were also divided into blocks of 52 representing 19-23 
flight hours each. Ten runs were made on each set of data 
and by using the mean damage value each damage calculation 
was multiplied by the respective number of blocks to produce 
the damage accrued for 1000 hours. The same procedure was 
followed for a Miner's Rule prediction. The results are 
shov/n in Table XV. 

Increasing the number of blocks from 26 to 52 resulted 
in an increase in the damage life for 1000 hours of 17-9 
percent v.'hereas the Minor's damage prediction sliowed a 
negligible Increase of .396 percent. Tlie jx’ason for- tills 
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TABLE XV 



A-7B Damage Accumulated for 1000 Plight Hours 
'Based on Accelerometer Readings of 146, 

64 j l8, and 1 for 52 Blocks 



Method 


Nc 


1 . of Blocks 


Damage /Block 


Total Damage 


ESDU 




26 


.0008126 


.021128 


ESDU 




52 


.000479 


.024908 


Miner ' s 




26 


.00405351 


.1053913 


Miner ' s 




52 


.0020348 


.1058096 


Increase is 


due 


to the fact 


that everytime a 


block begins. 


the material 


is 


essentially 


free' of residual 


stresses, and 



thus no favorable mean stress correction occurs until a high 
g load causes plastic deformation. The Miner's approximation 
is not concerned with the build-up of residual stresses, and 
thus should show no significant change in damage prediction 
as the number of blocks increases. 

It was assumed that after every flight the favorable 
residual stresses were relieved by the damaging effect of 
landings and the taxiing phase . It should be noted that the 
taxiing phase will only have a damaging effect on fatigue 
critical points outboard of the main mounts. In addition it 
was also assumed that the damage prediction increased linearly 
as the number of blocks Increased. Based on these assumptions, 
a damage calculation for the 1000 flight hours was divided 
into 500 flights that would represent an average of 2 hours 
between actual residual stress relieving. The damage 
calculated based on a .linear relationship resulted in a 



damage prediction of .09 or an expected life of 11,000 hours. 
This was not considered totally satisfactory since NADC's 
predictions are producing life estimates [Ref. in excess 
of 22,000 hours (excluding the factor of safety). 

B. EFFECTS OF STRESS RELIEF 

It was decided to Investigate the relationship between 
increasing flight blocks and damage prediction.. The A- 7 B 
previously used for a one-block flight was divided into 
flight blocks of one, two, three, four, six, and twelve. 

The damage prediction was then calculated for the entire 
31.5 hours of flight time. The results, including Miner’s 
prediction, are indicated in Table XVI. Ten runs v^^ith each 
set of data were run with the mean value of the damage for 
the ten runs tabulated. 



TABLE XVI 



A- 7 B Damage Predictions with Increasing Flight 


Blocks 


No. of Blocks 


Damage/Flight 

ESDU 


Total 
Miner ’ s 


Damage 

ESDU 


1 


. 0017 ^^ 0^4 ^4 


.0117782^ 


. 0017 ^ 10 ^ 


2 


.00106697 


.01215122 


.0021339 


3 


.00070075 


.0019376 


.00210225 




. 000588 ^ 4^4 


.01185008 


.00235376 


6 


.000^3121 


.0123125 


.00025873 


12 


.00023218^4 


.01176695 


.0027862 



Except for the decrease in prediction from a block of 
three flights the damage Is clearly Increasing as should be 
expected. The predictions from Miner’s Rule are relatively 
stable. Although the Increase In damage prediction Is not 
linearly related to an Increase In flight blocks^ 
the relationship Indicated a linear Increase with the 
loga.rlthm of the number of flight blocks. Figure 17 Indi- 
cates the trend. Based on these results an aircraft damage 
prediction could be determined by using a one flight block 
run and an additional block run that would be compatible 
with Input data. Then, assuming a linear relationship between 
the damage and the logarithm of flight block, the damage for 
any number of flights during that block could be determined. 
Based on this relationship the damage for the 500 flight 
block of the oscillograph data for 1000 flight hours produced 
a damage prediction of .0372 per 1000 hours or a life expec- 
tancy of 26,851 flight hours. This may seem a rather high 
estimate, but the loading profile the aircraft was designed 
for Is considerably higher than the actual loads experienced. 

It Is not suggested that the linear relationship of the 
logarithm of the number of flights Is clearly a linear 
relationship over the entire range. Extending predictions 
based on 26 or 52 flights for 5 OO flight blocks would not 
produce accurate results. Hovjever if the number of flights 
was close to the blocks used in obtaining two point predic- 
tions the linear relationship should produce reasonable 
results. In the case of the A-7B with 31-5 hours of 



accelerometer data, blocks of one and twelve flights should 
predict a -reasonable damage prediction for the aircraft if 
in fact 15 flights were involved during this reporting period 
which results in an average flight time per flight of 2.1 
hours . 

How much relief of residual stresses is caused by the 
landing and taxiing portion of the flight is presently 
unknown. Assuming that all the residual stress is relieved 
as a result of the ground phase is as unreasonable as assuming 
that none of the residual stress built up by a, high load is 
relieved. Somewhere in between these two extremes is the 
actual amount of stress relief that occurs. Until this 
information is known the assumption that total stress relief 
occurs after each flight is a highly conservative estimate 
that would produce damage predictions that have a built-in 



factor of safety. 



IX. SUMMARY AND RECOMMENDATIONS 



A need has been demonstrated for the implementation of 
a damage hypothesis in fatigue prediction that accounts for 
sequential effects. Predictions based on block loading 
programs have shown to be an optimistic prediction of the 
fatigue dam.age . Recording devices that have a direct rela- 
tionship to the stress experienced at a fatigue critical 
point such as time history strain gauges are highly desirable. 
Due to the nature of aircraft loadings, most fatigue damaging 
loads are initiated from the one-g level or slightly belox'J 
it. As a result of this type of loading, accelerometers were 
considered acceptable for indicating the number of loads that 
occurred above the preset levels. Unfortunately no informa- 
tion v;as provided concerning the order in vjhich the loads 
occurred. It was found that the occui*rence of a particular 
load in flight could not be sequentially associated vjlth the 
occurrence of any other load. As a result it was concluded 
that maneuver loads occur randomly. 

Based on accelerometer data, a flight profile was devel- 
oped that includes all fatigue damaging loads that a 
particular type aircraft experiences. Unfortunately the 
load levels now being recorded on the A-6A and A-7B aircraft 
are above the level necessary to obtain the flight profile. 

The need for recording load levels above the eight-g 
level in A-7B aircraft is appi’eclated due to the lilghly 



137 



damaging effect these loads Impose; however, analysis of 
flight data indicates that loads above the eight-g level 
occur on the average less than five times for every 1000 
hours of flight time. Reference 13 Indicates that loads 
with a peak of 8.05 g's produce a damage of .0001. Based 
on these data an A-7B which experienced only 8.05 g's could 
fly for 2,000,000 hours before fatigue failure. It is 
recommended that the accelerometer levels on A-7B aircraft 
be lowered to Include readings at the four-g level. 

Estimation of elght-g occurrences could be provided by 
allowing a reasonable number of the seven-g loads to occur 
above the eight-g level by means of some distribution 
similar to the equations provided in Table III. 

Randomizing the order of flight loads has a significant 
effect on the damage calculation provided by the ESDU method. 
During the life span of the aircraft, it is anticipated that 
the majority of the loading sequences generated would have 
actually been experienced by the aircraft. Varying the 
ratio of loads in each range had a less significant effect 
on the damage prediction, since the Miner's Rule approximation 
produced approximately similar results regardless of the 
number of loads a particular trial run produced. Neglecting 
a variation of load level ratios would preclude the possibility 
of one aircraft accidentally being assigned an extreme low 
ratio for each flight period. 

The flight profile generated from a small sample of 
aircraft cannot be considered as conclusive evidence that 
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each aircraft has a "fingerprint" flight profile below some 
high load level. Oscillograph data now compiled by type of 
aircraft should be reduced on an individual aircraft basis 
so that the profile hypothesis can be substantiated. 

The ESDU cumulative damage hypothesis provides a reason- 
able explanation of the overconservativism of the Miner's 
approximation. The ESDU method is suitable for computer 
analysis of damage prediction of a large fleet of aircraft 
due to the short run time required for calculation. 

Assuming that the residual stress imposed on the air- 
craft component due to a high load would remain indefinitely 
is unreasonable. Hov;ever the amount of stress relief caused 
by landing and taxiing loads is not available. It is recom- 
mended that the amount of stress relief imposed by the ground 
phase of the flight be given considerable attention to enable 
an increased accuracy of dam.age prediction by the ESDU method. 
Until information is known concerning the ground phase of 
the flight, each flight could be assumed to relieve completely 
the favorable residual stresses that developed during the 
flight . 
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APPENDIX A 
PROGRAM 5SDU 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

PROGRAM ESDU UTILIZES THE ESDU 
CUMULATIVE DAMAGE HYPOTHESIS FOR 
PREDICTING THE EXPENDED FATIGUE 
LIFE GF A-7B A/C FROM ACCELEROMETER 
READINGS AT THE 4t5t6, AND 7 G LEVEL 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



DIMENSION UPPEAK(5000) , DNPEAK ( 5000 ) , SEQUESdOOOU 
WRITE (6,&J 
READ (5,72 NUM 

DO 5 

READ (5,81 NACRFT, FOUR, FIVE , SIX, SEVEN 

AJUST=1 .0 

PREDIC=0eO 

IF (F0UR.GT«325.02 GO TO 1 
GO TO 2 

1 CALL FIX (FOUR, FIVE, SIX, SEVEN, AJUST) 

ASSURES MEMORY REQUIREMENTS SATISFIED 



C 

C 



SUBROUTINE TO GENERATE MAXIMUM G PEAKS 
A/C HAS EXPERIENCED 

2 CALL UPSPEC ( NPEAK , NACRFT , TOUR , F I VE , S I X , 
ISEVEMjUPPEAK, II I n 
IF (NPEAKoEQ.OJ GO TO 4 



GENERATE MINIMUM PEAKS 



CALL DNSPEC ( NPE AK , NACRFT , DNPEAK , FOUR , 1 1 1 I ) 
IF (NPEAK.EQ.O) GO TO 4 



GENERATE ODD INTEGER FOR 
RANDOM GENERATOR 



NODD= ( I FI X( FO UR ) +I F IX { FI VE ) + 1 F I X ( 5 1 X ) -i- 
IIFIXI SEVEN! i=!=I 1 1 I 
IF (NGDD.EQ.Ol GO TO 4 
IF (M0D{N0DD,2) cEQ.l) GO TO 3 
NODO=-NODD.-H 



LOADINGS PLACED IN RANDOM ORDER 

3 CALL MIXER ( IJPP EAK , DNPEAK , S EQUE S , NPE AK , NODD ] 

OBTAIN FATIGUE DAMAGE 

CALL ESDU (SEQUES ,NPEAK,PREDIC) 
PRED1C='5R,6DIC>:'AJUST 

4 WRITE (6,9! NACRFT , FOUR , F I V E , S I X , SEVEN , P RE D JC 

5 CONTINUE 

STOP 



6 FORMAT ('0'.' N AC R FT • , T 14 , ' FOUR ' , T29 , • F I VE » , T44 , ' S I X ' , 
3.T59, ' SEVEN" ,T73, ' PREOIC' , // ) 

7 FORMAT (I 10) 

8 FORMAT ( 110, 'IF i 0.0 I 

9 FORMAT ( ' 0‘ ,10 , T 1 5 , F4 , 0 ,T 29 , F4. 0 , T43 , FA . 0 ,T60,F4. 0, 

IT /i,ro. 7,//] 

END 
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cccccccccccccccccccccccccccccecccccccccccccc-ccccccccccccccc 

PROGRAM ESOU UTILIZES THE ESDU 
CUMULATIVE DAMAGE HYPOTHESIS FOR 
PREDICTING THE EXPENDED FATIGUE 
LIFE OF A-7B A/C FROM ACCELEROMETER 
READINGS AT THE 4,5,6, AND 7 G LEVEL 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



DIMENSION UPPEAK(5000), DNP EAK ( 5000 J , SEQUESdOOOU 
WRITE (6,6) 

READ (5,7) NUM 

00 5 

read (5,8) NACRFT, FOUR, FIVE , SIX, SEVEN 

AJUST=1.0 

PREDIC=0.0 

IF (FOUR. GT. 325.0) GO TO 1 
GO TO 2 

1 CALL FIX ( FOUR, FIVE, SIX, SEVENvAJUST) 

ASSURES MEMORY REQUIREMENTS SATISFIED 



SUBROUTINE TO GENERATE MAXIMUM G PEAKS 
A/C HAS EXPERIENCED 

2 CALL UPSPEC (NPE AK , NACRFT , FOUR , F I VE , SI X , 
1SEVEN,UPPEAK,IIII) 

IF (NPEAK.EQ.O) GO TO 4 

GENERATE MINIMUM PEAKS 

CALL DNSPEC ( NPE AK , NACRFT , DNP EAK , FOUR , 1 1 1 1 ) 

IF (NPEAK.EQ.O) GO TO 4 

GENERATE ODD INTEGER FOR 
RANDOM GENERATOR 

NGDD=( IF1X{FGUR)+IFIX(FIVE)+IFIX(SIX)+ 

IIFIX(SEVEN) )^I I I I 
IF (NCDD. EQ.O) GO TO 4 
IF (M0D(NQDD,2) . EQ.l) GO TO 3 
NGDD=N0DD+1 

LOADINGS PLACED IN RANDOM ORDER 

3 CALL MIXER ( UPPEAK, DNPEAK, SEQUES, NPEAK, NODD ) 

OBTAIN FATIGUE DAMAGE 

CALL ESDU (SEQUES, NPEAK, PREDIC) 

PREDIC = PR£DIC=^AJUST 

4 WRITE (6,9) NACRFT, FOUR,FIVE,SIX, SEVEN, PREDIC 

5 CONTINUE 
C 

STOP 

C 

6 FORMAT (*0*,‘ N AC RFT ' , T 14 , * FOUR ' , T29 , * F I VE • , T44 , ' S I X * , 
1T59,'SEVEN‘,T73,' PREDIC ,// ) 

7 FORMAT (110) 

8 FORMAT (I10,4F10.0) 

9 FORMAT ( 'O' , I8,T15,F4.0,T29,F4.0,T43,F4.0,T60,F4.0, 
1T71,F8.7,//) 

END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C SUBROUTINE REDUCES DATA FOR COMPATIBLE 

C MEMORY STORAGE WITH DIMENSIONED ARRAYS 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



SUBROUTINE FIX ( FQ-JR, F I VE , S I X , S EVEN» AJUST ) 

X=325.0/F0UR 

F0UR=325.0 

FIVE=X'-^FI VE 

SIX=X=«=SIX 

SEVEN = X=i'SEVEN 

IF (AMDDIFIVE, 1. ) .GT.0,5) F I VE = FI VE+l . 0 
IF (AMODISIX, 1. ) .GT.0.5) SI X=SI X+1 .0 
IF (AMOD(SEVEN?ic ) *GT.0c5) SEVEN=SEVEN+1 .0 
FIVE=AINT( FIVE) 

SIX=AINT( SIX) 

SEVEN=AINT(SEVEN) 

AJUST=1 .0/X 

RETURN 

END 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C SUBROUTINE GENERATES MAX G PEAKS 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



SUBROUTINE UPSPEC ( NP E AK , NACRFT , FOUR , F I VE , S I X , S EV EN , 
lUPPEAKf II I I) 

DIMENSION UPPEAK(5000), NPEK(12) 



C 

C 

C 

C 

C 



USING NORMAL DISTRIBUTION CHOOSE NUMBER 
OF PEAKS IN EACH LEVEL 



DO 1 I=1j12 
I NPEK(I)--=0 



C 

C 

C NEED ODD INTEGER TO INITIATE CALL 

C 

NODD= I F IX ( FOUR ) -M F I X ( F I VE ) -i- 1 F I X ( S I X ) + 1 F I X ( S E VEN ) 
IF (NODDoEQ.O) GO TO 13 
NCDD=U'JACRFT/NODD)=MII I 
IF (M0D(N00D,2) .EQ.l) GO TO 2 
NGDD=NCD0+1 
2 IX=NODD 



C 

C FIRST RANGE 

C 

IF (IF! X( FOUR) .EQ.O) GO TO 12 
C 

C NORMAL OISTIRBUTION RANDOM NUMBER 

C GENERATOR 

C 

3 CALL GAUSS ( I X , . 9764, 4 . 01 , V ) 

IF (V.GT.5»62) GO TO 3 
IF (V.LT. 2.399) GO TO 3 
V=FOUR*V 
NPEKtl)=IFIX(V) 

C 

C ROUND OFF TO NEAREST WHOLE NUMBER 
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c 

IF (AMDD( V,l. ) .GT.0.5) NPEK ( 1 ) =NPEK ( 1 J + 1 
SECOND RANGE 

4 CALL GAUSS ( I X, . 372 , 1 . 92i V) 

IF (V.GT.2.533) GO TO 4 
IF (V.LT. 1.306) GO TO 4 
V=F0UR*V 
NPEK(2) = IFIX( V) 

IF (AM00(V,1. ) .GT.0.5) NPEK ( 2 ) =NPEK( 2 H- 1 
THIRD RANGE 

5 CALL GAUSS ( I X , . 5058, 1 . 65 » V ) - 
IF (V.GT.2.484) GO TO 5 
IF (V.LT. .8154) GO TO 5 
V=F0UR^'V 
NPEK(3)=IFIX( V) 

IF (AM0D(V,1. ) .GT.0.5) NPEK (3 ) =NP EK { 3 ) H-1 
FOURTH RANGE 

6 CALL GAUSS ( I X, . 314, . 947, V ) 

IF (V.GT. 1,465) GO TO 6 
IF (V.LT. .4289) GO TO 6 
V=FGUR*V 
NPEK(4)=IFIX(V) 

IF (AMOD( V, 1. ) .GT .0.5 ) NPEK ( 4 ) = NP EK ( 4 ) -«- 1 
FIFTH RANGE 

7 CALL GAUSS ( I X, . 1 0 1 , . 6315 , V ) 

IF (V .GT. .79815) GO TO 7 
IF (V.LT. .46485) GO TO 7 
V=F0UR*V 
NPEK(5)=--IFIX( V) 

IF (AM0D( V,l. ) .GT.0.5) NP EK ( 5 ) =NPEK( 5 ) + 1 
SIXTH RANGE 

8 CALL GAUSS ( IX, .0282, .417, V) 

IF (V.GT. .46353) GO TO 8 
IF (V.LT. .37047) GO TO 8 

V = FOUR='-V 
NPEK(5) =IFIX(V) 

IF (AHDD(V,1.) .GT.0.5) NP EK ( 6 ) =NP EK ( 6 ) + 1 
SEVENTH RANGE 

9 CALL GAUSS ( I X , . 0 195, .324, V ) 

IF (V.GT. . 3562) GO TO 9 
IF (V.LT.. 2918) GO TO 9 
V=FQUR>;'V 
NPEK(7)=IFIX( V) 

IF (AMC)D(V, 1. ) .GT.0.5) NPEK ( 7 ) = NP EK( 7 ) + 1 
IF (FOUR.LT.l.i) NPEK(7)=1 

EIGHTH RANGE 

10 CALL GAUSS ( I X , . 0 1 1 1 , . 263 , V ) 

IF (V.GT. .281) GO TO 10 
IF (V .LT. ,2447) GO TO 10 

V = FOUR'r^V 
NPEK(8) =IFIX( V) 

IF (AMOD(V,l. ) .GT.0.5) NP EK C 8 ) = NP EK ( 8 ) + 1 
IF (FGUR.LT.2. 1 ) NPEK(B)=1 

NINTH RANGE 

11 CALL GAUSS ( I X , . 0 1 7 5 , . 1 505 , V ) 

IF (V.GT.. 1794) GO TO 11 
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IF (V.LT..1216) GO TO 11 

V=FOUR*V 

NPEK( 9) =I FIX( V) 

IF (AMOD( V,l. ) .GT.0.5) NFEK ( 9) =NPEK( 9 ) + 1 

12 NPEK(10)=IFIX(FIVE) 

NPEK(11J=IFIX(SIX) 

NPEK(12)=IFIX(S£V£N) 

FILL VECTOR WITH LOADINGS 

13 CALL HI LOAD ( NPEK , UPPEAK, NP EAK ) 

RETURN 

END 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE ASSIGNS REPRESENTATIVE 
G LOADS TO EACH RANGE 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



SUBROUTINE HILOAD ( NPEK jUPPEAK , NPEAK J 
DIMENSION NPEK(12), UPPEAK(5000) 
NPEAK=0 
C 

DO 1 1=1,12 

1 NPEAK=NPEK{ I )+NPEAK 
C 

IF (NPEAK*EQ.O) GO TO 27 

NDUMMY=0 

NT0P=0 

IF (NPEK( U.EQ.Q) GO TO 3 
NDUMMY=NTOP+l 
NTOP=NDUHMY+NPEK( 1 )-l 
C 

DC 2 I=NDUMMY,NTOP 

2 UPPEAK(I)=2.25 
C 

3 IF (NPEK( 2) . EQ. 0) GO TO 5 
NDUMMY=NPEK(1)-!-1 
NTOP=WDUMMY+NPEK( 2)“1 

C 

DO 4 I=NDUMMY,NTOP 

4 UPPEAKtl ) =2.7 
C 

5 IF (NPEK(3 J . EQ.O ) GO TO 7 
NDUMMY=NTOP+l 
NTOP=NDUMMY-H\!PEK( 3 )-l 

C 

DO 6 I=NDUHMY jNTOP 

6 UPPEAKI I 1=3.1 
C 

7 IF (NPEKC 4) . EQ. 0) GO TO 9 
NOUMMY = NTOP-H 
NT0P=NDUMMY+NPEK(4)-1 

C 

DO 8 I=NDUMMY,NTOP 

8 UPPEAKI I 1 =3.45 

c 

9 IF (NPEK( 5) . EO.OJ GO TO 11 
NDUMMY=NTGP+1 
NTOP=^iDUMMY^-rJPEK( 51-1 

C 

DO 10 l=NDUMHY,NTOP 
10 UPPEAKd 1=3.75 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



11 IF (NPEK(6) . EQ.O ) GO TO 13 
NDUMMY = NTOP+-l 
NTOP=NDUMMY+NPEK ( 6 )-l 

DO 12 I =NDUMMY,NTOP 

12 UPPEAKI I )=4.05 

13 IF (NPEK(7 I.EQ.O) GO TO 15 
NDUMMY=NTOP+l 
NTOP=NDUMMY+NPEK( 7)-l 

DO 14 I=NDUMMY,NTOP 

14 UPP£AKm=4.35 

15 IF (NPEKI 8 ) . EQ.O) GO TO 17 
NDUMMY=NTOP+l 
NTOP=NOUMMY+-NPEK { 8 ) "1 

DO 16 I=NDUMMY,NTOP 

16 UPPEAKd )=4.65 

17 IF (NPEK(9).EQ.O) GO TO 19 
NDUMMY=NTOP+l 
NTOP=NDUMMY+NPEK(9)-l 

DO 18 

18 UPPEAKI I )=4.9 

19 IF (NPEK( lOJ .EQ.O) GO TO 21 
NDUMMY---NTOP + 1 
NTOP=I\'DUf>1MY+NPEK( 10 )-l 

DC 20 I=NDUMMY,NTOP 

20 UPPEAKU )=5«75 

21 IF (NPEK( 11) cEQ.O ) GO TO 23 
NDUP.MY = NT0P+1 
NTOP=NDUMMY+NPEK( 11)-1 

DO 22 I=NDUMMY,NTOP 

22 UPPEAKd )=6. 6 

23 IF (NPEK( 12) .EQ.O) GO TO 25 
ivlDUM.MY = NTOP^l 
NTOP=NDUMMY+NPEK( 12)-1 

DO 24 I =NDUMMY,NTOP 

24 UPPEAKd)=7,8 
END 
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CCCCCCCCCCCCCCCCCCCCCCCCCCCCC.CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C SUBROUTINE GENERATES MINIMUM PEAKS 

c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 



SUBROUTINE ONSPEC (NPEAK» NACRFT , DNPEAK, FOUR , 1 1 1 1 ) 
DIMENSION DNPEAK(5000) , NPEK(5) 

C 

DG 1 1=1,5 

1 NPEK( I ) =0 
C 

IF (IFIX(FCUR).EQ.O) GO TO 7 
NCOD=(NACRFT/FOUR)*IIII 
IF (M00(N0DD,2) . EQ.l) GO TO 2 
NOOD=NODD-H 

2 IX=NODD 
C 

C FIRST RANGE 

C 

3 CALL GAUSS (I X , 1 . 93099 , 7. 1554, V ) 

IF (V.GT. 10.342) GO TO 3 

IF (V.LT. 3.9692) GO TO 3 

V=FDUR*V 

NPEK(1)=IFIX(V) 

IF (AM0D(V,1.) .GT..5) NPEK( 1 ) =NPEK ( 1 ) +1 
C 

C SECOND RANGE 

C 

4 CALL GAUSS ( IX, .8563, 1.6728 ,V ) 

IF (V.GT. 3. 0857) GO TO 4 

IF (V .LT, . 2599) GO TO 4 

NPEK(2)--=IFIX( V) 

IF (AM00(V,1.) .GT..5) NPEK{ 2)=NPEK( 2)+l 
C 

C THIRD RANGE 

C 

5 CALL GAUSS ( I X . . 1 8698 , . 3588 8, V ) 

IF (V.GT. .6674) GO TO 5 

IF (V .LT. .050363 ) GO TO 5 

V=FOUR«V 

NPEK(3)=IFIX(V) 

IF (AM0D(V,1.) .GT, .5) NPEK ( 3) = NP EK ( 3 ) -H 
C 

C FOURTH RANGE 

C 

6 CALL GAUSS ( I X, . 05951 10252 , V ) 

IF (V .GT. .200712 ) GO TO 6 

IF (V.LT. .00433) GO TO 6 
V = FDUR=i=V 
NPEK(4) =IFIX( V) 

IF (AM0D(V,1. ) .GT. .5) NPEKl 4)=NPEK(4) +1 
C 

C FIFTH RANGE 

C 

7 NPEK(5) =NPEAK-(NPEK(1 ) +NPEK ( 2 ) +NP EK( 3 ) +NPEK ( 4 ) ) 
CALL LOLOAD ( NPEK , DNPEAK) 

RETURN 

END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE ASSIGNS G LOADINGS TO 
MINIMUM PEAK RANGES 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCCCCCCCCCC 



SUBROUTINE LOLOAD ( NP EK , DNP EAK ) 
DIMENSION NPEK(5), DNPEAK(5000) 
NPEAK=0 
C 

DO 1 1=1,5 

1 NPEAK=NPEK( n+NPEAK 

c 

NDUMMY=0 

NT0P=0 

IF (NPEK( 1) . EQ.O) GO TO 3 
NDUMMY=NTOP+l 
NTOP=NDUMMY+NPEK( 1 )-l 
C 

DO 2 I=NDUMMY,NTOP 

2 DNPEAKI I ) = .538 
C 

3 IF (NPEK(2) .EQ.O) GO TO 5 
NDUMMY=NTOP+l 
NTOP=NDUMMY+NPEK(2 )-l 

C 

DO 4 I=NDUMMY,NTOP 

4 DNPEAKU )=.346 
C 

5 IF (NPEK(3) .EQ.O) GO TO 7 
NDUMMY=NTOP+l 
NTOP=NDUMMY+NPEK ( 3 )-l 

C 

DO 6 I=NDUMMY,NTOP 

6 DNPEAKd ) = .148 
C 

7 IF (NPEK(4) . EQ.O) GO TO 9 
NDUMMY = NTOP-M 
NT0P=NDUMMY-s-NPEK(4)-l 

C 

DC 8 I=NDUMMY,NTOP 

8 DNPEAKU )=--.054 
C 

9 NDUMMY=NTOP+l 
NTOP=NDUMMY-i-NPEK ( 5)-l 

C 

DO 10 I =NDUMMY,NTOP 

10 DNPEAK ( I )=1.0 
C 

IF (NPEAK.NE.NTOP) GO TO 11 
GO TO 12 

11 WRITE (6,13) 

NPEAK=0 

12 RETURN 
C 

END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C SUBROUTINE GENERATES UNIFORMLY 

C DISTRIBUTED RANDOM NUMBER 

C 

cccccrxcccccccccccccccccccccccccccccccccccccccccccccccccccc 



SUBROUTINE RANDU (iXtIY,YFL) ^ 

IF (lY) 1,2,2 

1 IY=IY+2147483647+1 

2 YFL=IY 

YFL=YFL«. 465661 3E-9 

RETURN 

END 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C SUBROUTINE GENERATES NORMALLY 

C DISTRIBUTED RANDOM NUMBER 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 



SUBROUTINE GAUSS (IX,S,AM,V) 

A = 0.0 
C 

DO 1 1=1,12 
C 

C UNIFORMLY DISTRIBUTED RANDOM 

C NUMBER GENERATOR 

C 

CALL RANDU ( iX, lY, YJ 
IX=IY 
1 A=A+Y 
C 

V=(A-6. 0) *S+AM 

RETURN 

END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C SUBROUTINE ARRANGES PEAKS IN 

C A RANDOM ORDER 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC.CCCCCCCCCCCCCCCCCCCCC 



SUBROUTINE MIXER ( UPP EAK, DNPE AK , S EQUES , NP E AK , NCDD ) 
DIMENSION UPPEAKI5000 ), DNP EAK ( BOOO ) , SEQUESUOOOU, 
1M(5000) 

NOTEVN=NODD 

CALL SHUFLE ( M , NP EAK, NOTE VN ) 

C 

C RANDOMIZE INDEX OF LOADS 

C 

c 

DO 1 I--ltNPEAK 

1 SEQUES( Z-I-1 J = DNPEAK(M( I) ) 

C 

NOTEVN=NODD+8 

CALL SHUFLE ( M » NP E AK, NOTEVN ) 

C 

DC 2 I=1,NPEAK 

2 SEQUES( 2':'I )=UPPEAKIM{ I) ) 

C 

SEQUESC 2=«'NPEAK + 1)=1.0 

RETURN 

END 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C SUBROUTINE SHUFFLES INDEXES CF M 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



SUBROUTINE SHUFLE ( M, NPEAK, NOT EVN ) 
DIMENSION MINPEAKI 
CALL OVFLOk 
C 

DO 1 1=1, NPEAK 
M( I ) = I 

1 CONTINUE 
C 

K=NPEAK 

2 CALL RANDOM (NOTEVN,U,l) 

I=K*U-5-l 

L=H( I ) 

M(I ) = M(K) 

M ( K ) = L 

IF (K.GT.l) GO TO 2 

L=M(NPEAKJ 

M(NPEAK I=M(2i 

M(2) = L 

RETURN 

END 
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ccccccccccccccccccccccccccccceccccccccccccccccccccccccccccc 

SUBROUTINE CALCULATES FATIGUE LIFE 
BY ESDU METHOD 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 



SUBROUTINE ESDU ( SEQUES ,NPEAK, PREDIC) 

DIMENSION SEQUES( 10001) 

CORECT=0.0 

TK=3.5 

B0T=3.6 

SLITLE=5.5 

FT=77.0 

FP=69.0 

SMELCA=10c0 

PREDIC=0.0 

NMCST = 2>i'NPEAK + l 

DO 1 I=1,NM0ST 

CONVERT TO STRESSES IN KS I 
1 SEQUESI I )=3.5143*SEQUES(I ) 

NMOST=NM0ST-l 

COMPUTE STRESS MEAN ACTUAL PROGRAM 



00 16 

0RRECT=0.0 

SM-( SEQUES' I + l) ^SEQUESf I » »/2*0 
SA=ABS{ SEQUES ( I U -SEQUES ( I ) )/2.0 

CONVERT TO ACTUAL STRESS 

SMEL=TK=«=SM 

SAEL=TK«SA 

SMAXEL=SMEL+SAEL 

SMINEL=SMEL-SAEL 

CHECKING FOR TYPE 4 

IF (SAEL.GT.FP) GO TO 2 
GO TO 3 

2 IF (SEQUES(n-l) .GT.SEQUES(I) j GO TO 4 
IF (SEQUESn+l) «LT.SEQUES{ I ) ) GO TO 6 

3 IF ( ( CORECT^-SMAXEL) vGTcFP) GO TO 4 
GO TO 5 

TYPE 3 

4 CQRECT=FP-SMAXEL 
GO TO 7 

5 IF ( ( CORECT + SMINEU «LT, {-l.=^^FP) ) GO TO 6 
GC TO 7 

TYPE 2 

6 CORECT=-l ,0-( SMINEL+FP) 

7 SMACP=SMEL+CCPECT 
SMAXEC = SMELCA^-SAEL 

DOWN LOADS NO CONTRIBUTION 
TO DAMAGE 

IF (SEQUESn-s-1) cLT.SEQUESn ) ) GO TO 16 
ADJUST FOR S N DATA 
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IF ( SMAXEC.GT.FP ) ORRECT=FP-SMAXEC 
SMACA=0RRECT4-SMELCA 
B0TCAL=1. 0-(SMACP/FT)**2 
TCPCAL=^ L . 0-{ SMACA/FT ) 

IF (SMACA. LT.0.0 ) TOPCAL= 1 . 0+ ( S MACA/ FT ) **2 
IF (SMACP.lt. 0.0) 80TCAL=1.0+(SMACP/FT)**2 
SAELCA= (T0PCAL/80TCAL )*SAEL 
SAELCA=SAELCA/TK 
IF (S AELCA.lt. BOT ) GO TO 16 



S N DATA 



CYCLES=9.8596-.5981-SAELCA+.0142 07«(SAELCA*’i‘2 ) 
CYCLES=10^=!'CYCLES 
PREOIC=1.0/CYCLES+PR£DIC 
16 CONTINUE 
C 

RETURN 

END 
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APPENDIX B 



Figures 

Three Different Loading Spectrums 

The Delaying Effect of Peak Loads on Crack 
Propagation 

Restricted Level-Crossing Method 

Differences Between Gust and Maneuver Spectra 

Gust Loading Spectrum and Maneuver Spectrum 

A-7B,^A-6A, A-^E, F~^B, and F-8E Exceedances 
Per 1000 Hours 

Ratio of 2-2.5 g Loads to Other g Levels 

Ratio of g Level Loads to 4-4.5 g Loads 

Ratio of Loads to 4 g Exceedances 

A-7B Exceedances Per 1000 Hours 

A-7B Ratio of Loads to 4-5 g Loads 

A-7B Proposed Ratio of Loads to 4-5 g Loads 

A-7B Ratio of Negative Loads to 4-5 g Loads 

A-7B Load Sequence 

A-7B Load Sequence 

A-7B Load Sequence 

A-7B Load Sequence 

Type 1 Load Cycles 

Type 2 Load Cycles 

Type 3 Load Cycles 

A-7B Damage Prediction vs. Flights 
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PIG. ia THREE DIEFEREKT LOADING SPECTRUMS 
APPLIED TO SPECIMENS^ 



19.2 

~2.9 

19.2 
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$c«t- CRACK LENGTH 
S (mm) 




7'l 



ALCLAD SHEET SPECIMENS, WIDTH 160 mm, THICKNESS 2 mm ( REF. 27). 






I 



O FIRST COUNTING CONDITION SATISFIED 
o SECOND »» tt tr 





PIG. 2 RESTRICTED LEVEL-CROSSIRG MSTHODo 
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PIG. 3 DIFFERENCES BETV/EEN GUST AND MANEUVER SPECTRA 
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PIGc 4 GUST LOADING SP3CTRUM (TOP) AND 
MNSUVER SPECTRUM (BOTTOM). 



\ 
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LOAD .( g * s ) 

PIGe 5 A-7B, A-6A,A-4E,P-4 b,AND P-8E 
EXCEEDANCES PER 1000 HOURS. 
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2.75 



4.75 



3*25 3«75 4«25 

MEAN OP LOAD LEVELS (g's) 

PIG. 6 RATIO OP 2-2.5 g LOADS TO OTHER LEVELS. 
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MEAN OF LOAD LEVELS (g's) 

PIG» 7 RATIO OP LOAD TO 4-^-. 5 e‘s FOR FIVE A/C. 
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PIG. 8 RATIO OP LOADS TO 4 g EXCEEDANCES . 
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*2 4 6 

LOAD (g*s) 

FIG. 9 A-7B EXCE5DANC3S PER 1000 HOURS. 
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LOAD (g's) 

PIG. 10 A-7B RATIO OP LOADS TO 4-5 g LOADS. 
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Ratio of Loads to 4 g Loads 



VO 



ir» 






OJ 




1234 

LOAD (g's) 



•PIG.n A-'TB PR0P0S3D RATIO OP LOADS TO 
• 4-5 g LOADS. 




10 AD (g's) 



fig, 12 A-7B PROPOSED RATIO 0? NEGATIVE 
LOADS TO 4-5 g LOADS. 
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TIME 



,I'IG.13 A-7B LOAD SEQUENCE BASED ON ACCELEROMETEH 
•READINGS OP 3,5,3, AND 




TIME 

.PIG.14 A-7B LOAD SEQUENCE BASED ON ACCELEROMETER 
READINGS OP 3,5,3, AND 4. 
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TIME 

PIG* 15 A-7B LOAD SEQUENCE BASED ON ACCELEROMETER 
READINGS OP 5,5,3, AND 4* 
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TIKE 

PIG,16 A»7B LOAD SEQUENCE BASED ON ACCELEROMETER 
READINGS 0? 3? 5? 3, AND 4. 
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fig, 17 TYPE 1 LOAD CYCLES. SMAXSL 
LESS THAN PP. 
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PIG. 18 TYPE 2 LOAD CYCLES. SMIKEL 
LESS THAN -PP. 
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FIG0I9 TYPE 3 LOAD CYCLES. SMAXEL 
GREATER THAN PP. 
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Fatigue Damage times 1000 

2*1 2*3 . 2,5 . 2.7 




2 4-8 16 

Flights 



PIG. 20 A-7B DAMAGE PREDICTION VS PLIGHTS 
(31.5 HRS.). 
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